A plastic flow theory for amorphous materials 
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Starting from known kinematic picture for plasticity, we derive a set of dynamical equations 
describing plastic flow in a Lagrangian formulation. Our derivation is a natural and a straightforward 
extension of simple fluids, elastic and viscous solids theories. These equations contain the Maxwell 
model as a special limit. We discuss some results of plasticity which can be described by this set of 
equations. We exploit the model equations for the simple examples: straining of a slab and a rod. 
We find that necking manifests always itself (not as a result of instability), except if the very special 
constant-velocity stretching process is imposed. 
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Plastic materials exhibit several features which are not 
present in the usual liquids or solids. Their dynamics con- 
sist in a nontrivial mixture of liquid-like and solid- like be- 
haviors. Understanding plasticity in metal industry and, 
in general, in technology, is of a paramount importance. 
Nonetheless, to date no universal dynamical equations 
describing plastic materials like Navier-Stokes equations 
for fluids, and Lame equations for elastic solids, are avail- 
able. Under strain, plastic material may exhibit elastic- 
like behaviors, yield stress, flowing behaviors, nonlinear 
engineering strain-stress relation, and so on [lj. A major 
goal in material science is the description of these phe- 
nomena in terms of dynamical evolution equations for 
relevant variables, namely the velocity, stress, and the 
analogues of strain. 

There has been important contributions to the theory 
of plasticity, especially for crystalline materials in terms 
of dislocations 0, [H, 0]. However, there is no need to 
evoke dislocations (if ever this notion has a meaning) for 
amorphous materials, and thus the question arises of how 
a corresponding theory can be build at the continuum 
level. This question has known recently an upsurge of 
interest [f| . An essential issue when addressing the ques- 
tion of plasticity is the distinction between crystalline 
solids and amorphous materials. Elastically deformed 
monocrystals are in a metastable state. Their plastic 
flow takes place only upon creation of dislocations, and 
is thus a nonlinear process. Conversely, plastic flow of 
amorphous materials should occur, in principle, linearly 
with respect to the applied stress. In crystals, an ad- 
ditional field, namely dislocation density, is introduced 
which couples to the elastic as well as to plastic distor- 
sions [1, IH • It is proposed here that one can derive a plas- 
tic continuum theory for amorphous materials, without 
evoking neither dislocation density, nor an internal vari- 
able that is distinct from the variables describing usual 



kinematics of plasticity. The evolution equations can be 
written in a closed form in terms of the elastic and plas- 
tic distorsions only. The concept of distortion used 
here was introduced by Kroener and Rieder in 1956 0], 
and will constitute our basic definition of the plastic flow 
variable. 

If Uk denotes the k th component of the displacement 
field, then the total distorsion tensor reads diUk- For a 
purely elastic solid, and in the small deformation regime, 
the symmetrized part — (diU^ + dkUi)/2 is the strain 
tensor. If plastic flow is involved, the total distorsion 
tensor d^k is a sum of the plastic flow contribution w?l, 



and the elastic part Wik (see Q) 



d l u k 



= w ik 



(1) 



The symmetrical part of the elastic distortion defines the 
strain tensor 
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u ik = -^{wih + w k i). 



(2) 



The constrain w\l = is usually supposed. Then from 



Eqs. (|ll2|) it follows u u = <9jMj. 

Usual elasticity can be presented in a Lagrangian man- 
ner by introducing an energy and a dissipative function. 
The elastic energy of a solid is given by 



A 2 
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(3) 



here A, \x are Lame coefficients, and the kinetic energy 
is pv 2 /2, where p is the density of the material. The 
dissipation function reads for a viscous solid 
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(4) 
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In the presence of a plastic flow one must introduce a new 
additional dissipative part related to plasticity. The sys- 
tem is described by three independent tensors, namely 
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Uik, and w\ k which are the symmetric (+), and anti- 
symmetric (— ) parts of the plastic distortion tensor, re- 
spectively. From the basic kinematic relation ([1]) and the 
strain tensor definition © it follows 

w fk + = ^(diUk + d k Ui) - u lk . (5) 

We expect the dissipation to consist of a quadratic form 
of these quantities, that we write as 

2«4 + +l(^) 2 +7(^j J (6) 

This is the dissipation corresponding to plastic flow. Sta- 
bility criteria for the dissipation function enforce 77 s > 0, 
£ > 0, 7777 s > a 2 , 7 > 0. Note that there is only one di- 
latational viscosity constant £. All other constants are 
related with shear motions. As we will see below in the 
liquid limit the constant 77 is a usual hydrodynamic vis- 
cosity. 

The strategy now consists in performing variations of 
the total Lagrangian with respect to the independent 
variables. Variation with respect to <5u, with 8wf k = 0, 
yields, upon using, 25ui k — diSuk + d k Sui, the momen- 
tum conservation law 

pili = d k (Tik (7) 

with the stress tensor Oi k consisting of a sum of the usual 
elastic part as well as the dissipative part with the usual 
solid viscosity terms, and an additional plastic term 

o~ik = ^uuSik + 2fiu ik + 

+ 2r] s u lk +(c~^V s ^j Sikuu + 2aw?+ . (8) 

Variation with respect to 5w\ l k + , with Su = 0, and 
fiwi k ~ = (in that case 5ui k = ~Sw^ k + ) provides us with 

Oik = a ^2u ik - ^SikUii^j + 2rjwf k + , (9) 

here di k is the traceless part of a ik . 

Finally, variation with respect to Swf k ~, with <5u = 0, 
and 5w^l + = (in that case Suik = 0) leads to uiffT = 
(a direct consequence of the absence of dissipation for 
rigid rotation; note also that energy does not depend on 
that mode). 

An important remark is in order. Differentiating (fl"j) 
with respect to time one obtains 

dkVi + d k Vt = -2jf k + 2u ik (10) 

where jfl = — w\ l k + , is the plastic current. This equa- 
tion (see also 0), apart from the (conventional) minus 
sign in front of jf k , bears resemblance with Eq.(3) of 
Ref. There is, however, a fundamental difference. In- 
deed, Eq.(3) of Ref. @ uses the kinematic condition (jTQJ) , 



plus Hooke's law, where Uik is assumed to be related to 
the stress tensor by 

_ Oik P_ 

Ulk ~ 2n 2R- 

where p = —Ukk/2, is the pressure, and K and /x are the 
compressibility and the shear modulus (note that a 2D 
geometry is assumed in Ref. [5|). In the present study 
we do not postulate a priori a Hooke's relation, since 
both elastic and plastic contributions are embedded to- 
gether within the total distortion tensor diUk- The rela- 
tion between Uik and o~ik follows here as a consequence 
of the Lagrangian formulation, and the relationship be- 
tween these two quantities is provided by ([8]) (showing 
that a measure of the stress is a combination of elastic 
and plastic deformations). 

It is possible to express the plastic distortion tensor 
in terms of other quantities. From Eqs.([8][9]) we may 
express o~ik in terms Uj& and its time derivative. It is 
convenient to split the stress tensor cr^ into a traceless 
and a pressure-like term (actually the trace of o~ik)'- 

0~ik = ^ikCll + CTifc (11) 

The trace has a usual elastic (including the dissipative 
part) form 

ou = (3A + 2/j)u« + 3C«h. (12) 

The traceless parts of the stress tensor are connected with 
each other and with spatial gradients of the velocity by 
the following two relations 

(77 - a)a ik = 2r\i±u ik + 2(7777* - a 2 )u ik , (13) 

a lk + 2(77 - a)u ik = 77 (diV k + d k v t - ^^div^J . (14) 

The first relation is obtained by expressing ujf k from 
Eq. ([9]) and inserting the resulting relation into (JSJ. The 
second one follows from Eq. © by using relation ([5]). The 
set of Eqs. (|7I1HI14[) defines space-time evolution of dis- 
placement vector u, strain Uik and the stress tensor o~ik- 
This constitutes a complete set of equations for the three 
(vectorial and tensorial ) quantities u, Uik and (we 
could, of course, alternatively use other quantities like 
w il + )- Note that Eq. (Tj"4"| has some similarity with the 
Maxwell model, used to describe plasticity with a yield 
stress in some models [7]. There is an important differ- 
ence, however. Instead of the &ik on the l.h.s. we have 
Uik- Again this consistently follows from the Lagrangian 
formulation. 

The Maxwell model of liquids with high viscosity can 
be obtained from our equations only if one consider the 
incompressible limit and set a = 77 s = £ = 0. This leads, 
from Eq. (fT3"|) , to cr,fe = 2fiUik- Eq. (fT4"|) reduces to the 
well known Maxwell form 

Oik + A"7~Vj fc = ^{diU k + d k Ui)- (15) 
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Let us present few examples where we could obtain 
an exact solution of the plastic dynamics. Consider 
an induced oscillatory motion in the material. We as- 
sume a semi-infinite medium bounded by a planar sur- 
face which undergoes oscillations in its own xz— plane: 
u x (t,0) = ucosujt. In this case the nonzero components 

and o~ xy . Then 



of the fields are u x , u xy , w xy , 
the set of Eqs. (j7ll3H4[) reads 



w 



pl+ 

xy 



PU X = dy(J X y = d,. 



2npu xy + 2 (7777 s 



L xy _ 



77 — a 



pu xy + (77 - 2a + r) )u xy = — - — d„v. 



y u x • 



(16) 



(17) 



It is then found that each field is a linear combination 
of the complex modes cx exp(— iuit + iny), where re is de- 
fined by 



(re'- 



■IK 



ll\2 



.uop 



1 - 



iu)(n — a)* 



77/j — iLo[j]r\ s — a 2 ) 



(18) 



For example, the displacement field is 
u x = u cos(ut — n'y) exp(— n"y). The low frequency 
limit recovers a known Stokes result for a shear viscous 
mode in liquids (see §24 in fgj]). Elastic solid behavior 
(an emission of shear sound) corresponds to the limit of 
high plastic viscosity 77 — > 00. 

We would like to point out some results that can be 
captured analytically in some special limit. The long 
time behavior of a slab under tension, is expected to be 
dominated by plastic flow. Ultimately, the plastic flow 
should look-like a hydrodynamical flow. Let us concen- 
trate on this limit. Consider a plate (or a rod) of a plastic 
material with free surfaces (FigJT]). This is a similar ge- 
ometry to that treated in Ref. Q ■ The plate is stretched 
along the x direction. For a flat geometry we have ob- 
tained an exact solution with the plate thickness h(t) 
that depends only on the t variable. This type of solu- 
tion exists only in the case where the stretching occurs 
at a given constant velocity. Let us first motivate the 
solution on the basis of symmetries. Because of the axial 
symmetry with respect to the y axis at x = 0, v x must be 
zero on that line. For constant h there is a simple solu- 
tion that fulfills that symmetry, v x = cx, where c(t) is for 
the moment an arbitrary function of time. From incom- 
pressibility condition we have v y — —cy + g{x,t), where 
g is a priori an arbitrary function of x and t. Symmetry 
with respect the middle line y = 0, enforces g = 0. 

We straightforwardly obtain from the Navier-Stokes 
equation the pressure field 



P 



,>[< +c 2 )y+p( C - C 2 )^ + /(<), 



where f(t) is a function of time to be determined below. 
At the free surface the normal component of the stress 
(the tangential vanishes automatically) must vanish. 
This is easily computed from the above result by using 



the definition a v 



—p + 2ndyV y — —p — 2r]c. Imposing 
Oyy = on the free surfaces at y = ±h, at any x, we 
obtain c + c 2 = 0. This provides us with c = (t — to) -1 , 
where to is a constant of integration. It is convenient 
to measure the time from the moment ioi so we will set 
to = 0. For the length of the strip L, one has L = cL, so 



L(t) 



FIG. 1: A sketch of the geometry under consideration. 

that L — st, with constant velocity s. Due to mass conser- 
vation one obtains that h = 7/t, where the total volume 
is defined as 7s. The solution corresponds to stretching 
plastic flow for t > 0, and s,j > 0, while for t < 0, and 
s, 7 < 0, it represents a contracting flow. 

Reporting this solution into the above- 
mentioned boundary conditions fixes the function 
f(t) = j 07 2 t -4 — 2 r qt~ 1 . The xx-component of the stress 
tensor is o~ xx = —p + 2r\d x v x = —p + 277c. The total force 
of stretching is defined as 



F, 



2d 



1 4 7 2 /i 8 V h 

^ dya xx = + — • (19) 



The first part is inertial, while the second one is vis- 
cous. The viscous part dominates at long time such that 
t 3> 7 2 / 3 z^~ 1 / 3 , where v = r\j p is the kinematic viscosity. 

It is a simple matter to extend the calculation to a 
cylindrical geometry (stretching of a rod) following pre- 
cisely the same line. We only give the results: The ve- 
locity is given by v z — cz, v r — —cr/2, with c(t) = t~ x . 
The pressure reads 



3p (l 



8t 2 \ t 



We obtain for the length and the radius L = st, 
a = j/yt. The volume is given by ira 2 L — irsj 2 . The to- 
tal axial force in the rod stretching problem is 



3d 



Itt [ 
Jo 



rdra zz = 



16 t 3 



S- 



(20) 



where S — ira 2 is a rod cross section area. 

The above solutions exists only for a constant velocity 
stretching. The question thus naturally arises of what 
happens if an other process is imposed. This is what we 
would like to investigate now. Following Ref.Q, if the 
lateral boundary of the plate moves at a pre-determined 
strain rate L/L = £1 = const., then our result shows 
that a homogeneous thinning of the strip is not possible. 
Thus a modulated strip prevails. This is a precursor of 
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the necking problem. Thus necking appears here as nat- 
ural phenomenon due to material flowing [i| whenever 
the stretching is not performed at a constant speed. The 
necking is not related with an instability Q, but rather 
the fact that a homogeneous thinning does not exist (ex- 
cept if L(t) oc t) in plastic dynamics. 

Let us investigate the stability of the homogeneous so- 
lutions. We consider deviations with large wavelengths 
as compared with the layer thickness. By analogy with 
the theory of shallow water (see [1| §108) we derive effec- 
tive hydrodynamic equations in terms of the thickness h, 
and velocity v along x. Mass conservation yields 

h + d x (hv) = 0. (21) 

Note, that when the viscous term dominates, one can 
rewrite the total force (|19p as follows i 7 ^ = —8ndth. In 
a general case of inhomogeneous shape one should have 
Fid = —8rj(h + vd x h). Consequently the momentum con- 
servation law in this purely viscous limit has the form 

h{v + vd x v)=-8vd x {h + vd x h). (22) 

From Eqs. (|21l22j) one obtains upon lin- 
earization about the 2D homogeneous so- 
lution: tSh + xd x 8h + ^fd x Sv + 5h = 0, and 
^(tSv + xd x 5v + Sv) = —8vtd x (t8h + xd x 5h). Intro- 
ducing the new coordinate x — x/t one arrives at: 
t 2 5h + 7<5w' + tSh = 0, and j(tSv + Sv) = -8vt5h'. 
Taking the Fourier transform oc exp(igi), we find for a 
given mode with wavenumber q of one of the fields, say 
Sh : t 2 Sh + A(t + T)5h + 2Sh = 0, here r = 2vq 2 . This 
equation has a first integral t 2 Sh + 2t5h + ArSh = C, 
where C is a constant. One obtains finally 




where t* is a second constant of integration. At 
large time, i>r, the deviation amplitude is 
5h = Ct _1 {l - 4-rt- 1 ln(t/et*)}. It decays mostly 
as t^ 1 like the strip thickness h. Consequently, one 
may say that the solution is marginally stable for large 
wavelength fluctuations. 

The 3D problem of the rod stretching with small in- 
homogeneity can be formulated in the same manner 
fEqs. (|21l22[) h The results are identical to the 2D ones, 
if one makes the substitutions h — > S, and 8v — > "Sv. 

We have solved Eqs. pT1l - ([22|) numerically in the fully 
nonlinear regime. We give here the major results: (i) If 
one imposes a constant stretching we find that the ul- 
timate stage is a homogeneous thickness that decreases 
in time as 1/t, in agreement with our analytical results, 
(ii) Starting from a small perturbation (of sine type), 
we observe marginal stability, (iii) Most importantly, we 
have found that if the stretching velocity is not constant 
the ultimate stage is a modulated thickness, of necking 
type. In fact, Fig[T] represents the result of our numeri- 
cal solution, that exhibits necking in the case of initially 
flat plate and L/L = const, (same stretching law as in 
Ref. Q). This behavior is found for various initial condi- 
tions, and (non constant) stretching laws. Thus necking 
seems to be a robust feature, which takes place whenever 
the stretching is non constant. 

It should be mentioned that here our plastic equations 
have been written by disregarding the so-called objective 
derivative (we have used ordinary derivatives) for ten- 
sors. One alterative in order to confer an objective form 
to these equations, is to replace the time derivative of 
tensors by the so-called co- rotational derivative [l(|, as 
is done in We shall report on full numerics of our 
completed set of equations (|7I1 ltil4|) in the future. 
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